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-  SUMMARY 
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The  four-axis,  elevation  over  ( azimuth  over  pitch  over 
roll,  director  mount  is  analysed  ^to  determine  the  torques 
generated  by  mass  unbalances  about  the  motor  driven  axes. 
The  effect  of  these  on  the  motion  of  the  stabilised 
platform,  and  the  effect  of  this  motion,  as  well  as  that 
of  mechanical  tolerances,  on  the  accuracy  of  direction  is 
examined.  The  analysis  is  restricted  to  the  case  of  the 
fully  erect  system,  firmly  connected  to  a  stationary, 
level  surface,  but  the  method  of  approach  could  be 
extended  to  more  complex  situations.  The  results  provide 
the  basis  for  the  design  of  small,  lightweight,  high 
performance  directors.  /  — 
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1.  INTRODUCTION 

A  necessary  component  of  many  weapon  systems  is  the  stabilised  mechanical 
mount,  used  for  directing  the  boresight  of  sensors  and  designators  (optical, 
infrared,  radar,  etc)  and  the  launching  rails  of  weapons,  (guided  missiles, 
guns,  etc).  It  must  be  able  to  point  the  boresight  accurately  in  any 
direction  (within  a  predetermined  sector  of  three  dimensional  space)  despite 
the  inevitable  changes  in  attitude  adopted  by  the  carrying  vehicle.  Two  axes 
of  motion  between  vehicle  and  boresight  is  the  minimum  number  required,  but 
practical  considerations  (such  as  the  difficulty  of  providing  unlimited  motion 
about  the  axes,  and  the  desire  to  reduce  the  inertia  of  rapidly  slewing 
structures)  often  result  in  the  use  of  three  or  four  axes.  Perhaps  the  most 
commonly  used  arrangement  is  that  of  the  elevation-over-azimuth-over-pitch- 
over-roll,  four-axis  mount.  This  report  considers  only  that  type. 

In  many  applications  (particularly  for  sensors  and  designators)  it  is  of  great 
importance  that  the  mount  be  small  in  size  and  light  in  weight,  but 
nevertheless  exhibit  high  pointing  accuracy  and  rapid  slewing  of  the 
boresight.  These  normally  conflicting  requirements  can  only  be  achieved  if 
the  design  uses  accurately  fitted,  light  but  strong,  components,  assembled 
into  well-balanced  structures,  driven  by  servo  motors  of  small  size  and  weight 
and  capable  of  only  modest  torque.  If  the  requirements  for  fitting  accuracy 
and  static  and  dynamic  balances  are  tighter  than  those  achievable  under  normal 
manufacturing  methods,  it  may  be  necessary  to  provide  for  ad justment -on-test 
of  these  parameters.  In  any  case  it  is  essential  that  the  necessary 
conditions  for  balance  and  accuracy  are  clearly  understood.  This  report  is 
concerned  with  an  analysis  of  these  conditions. 

Emphasis  is  given  in  this  report  to  the  particular  configuration  of  two 
nominally  identical,  separately  elevatable,  sensors  mounted  symmetrically  on  a 
single  structure,  which  rotates  at  a  constant,  significant  angular  rate  about 
the  (vertical)  azimuth  axis.  Analysis  is  restricted  to  the  case  of  operation 
on  a  stationary  level  vehicle,  with  small  deviations  from  the  vertical  of  the 
azimuth  axis.  This  restriction  results  in  mainly  linear  equations  of  motion 
which  are  directly  solvable,  but  still  gives  a  good  in-sight  into  the  problems 
of  balance  and  accuracy.  The  approach  could  be  extended,  however,  to  cover 
other  cases,  such  as  the  erection  process,  or  operation  on  a  moving  vehicle. 

The  analysis  begins  in  Section  2,  where  the  ideal  four-axis  mount  is 
described,  followed  by  the  definition  of  a  sequence  of  five  Cartesian 
coordinate  systems  attached  one  to  each  of  the  five  separate  solid  structures 
interconnected  by  the  four  axes  of  roll,  pitch,  azimuth  and  elevation.  Matrix 
equations  are  derived  relating  the  coordinates  of  a  single  point  referred  to 
each  of  the  five  coordinate  systems. 

In  Section  3  we  derive  general  relations  for  the  torque  components,  referred 
to  one  coordinate,  system,  caused  by  the  reaction  force  of  a  point  mass,  fixed 
in  a  second  coordinate  system  (neither  system  being  inertial),  and  subjected 
to  gravitational  and  kinematic  accelerations. 

In  Section  4  we  define  ten  mass  distribution  parameters  (moments)  for  each  of 
the  elevating,  rotating  and  combined  elevating  and  rotating  assemblies.  We 
calculate  torque  components  about  the  elevating  and  rotating  coordinate 
systems  due  to  each  of  these  parameters  under  varying  angles  of  roll,  pitch, 
azimuth  and  elevation.  These  torques  are  transformed  to  the  motor  driven  axes 
of  roll,  pitch,  azimuth  and  elevation,  and  the  results  are  presented  in 
tabular  form. 

Section  5  is  concerned  with  the  conditions  of  mass  distribution  for  static  and 
dynamic  balance  of  the  structure  (no  torques  about  the  motor  driven  axes). 
The  azimuth  rate  is  assumed  constant,  and  the  elevation  rate,  roll  angle  and 


WSRL-0355 -TR 


2 


pitch  angle  are  assumed  to  be  zero.  Ten  separate  parameters  must  be  zeroed. 

Sections  6,  7  and  8  are  concerned  with  the  behaviour  of  the  roll  and  pitch 
error  correction  loops,  under  the  assumption  of  constant  azimuth  rate,  zero 
elevation  rate  and  small  pitch  and  roll  errors.  The  servo  motor  torque 
demands  are  assumed  to  be  proportional  to  error  and  error  rate,  with  a  small 
unbalance  permitted  between  the  two  loops. 

It  is  shown  in  Section  6  that  a  uniform  spiral  motion  of  the  azimuth  axis 
about  the  vertical  position  generates  a  reaction  torque  (via  the  servo  motors 
and  the  inertia  parameters)  which  consists  of  three  spiral  components.  The 
main  component  is  at  the  same  angular  rate,  and  in  the  same  direction,  as  the 
motion.  One  minor  component  (at  much  smaller  amplitude)  is  at  the  same 
angular  rate,  but  in  opposite  direction,  and  is  proportional  to  the  unbalance 
between  the  servo  loop  gains.  The  second  minor  component  (also  at  a  much 
smaller  amplitude)  has  an  angular  rate  and  direction  depending  on  those  of 
both  the  spiral  motion  and  of  the  azimuth  scanning  rotation,  and  is 
proportional  to  the  magnitude  of  one  unbalanced  mass  distribution  parameter. 
The  magnitude  of  all  three  components  depends  on  the  angular  rate  of  the 
spiral  motion;  and  that  of  the  main  and  the  second  minor  component  depends  on 
the  relative  direction  of  the  spiral  and  scanning  motions. 

The  transient  response  of  the  roll  and  pitch  system  is  motion  in  the  absence 
of  any  externally  applied  forcing  torque;  and  hence  motion  producing  no 
reaction  torque.  This  is  calculated  in  Section  7  by  an  iterative  process. 
The  spiral  motion  giving  no  main  component  of  reaction  torque  is  determined. 
This  solution  is  then  perturbed  to  correct  it  for  the  minor  torque  components. 
An  infinite  series  of  spiral  motions  is  produced  by  this  method,  but, 
providing  the  unbalanced  mass  and  servo  gains  are  small,  the  amplitudes 
converge  rapidly  towards  zero. 

The  residual  (non-reaction)  torques  in  the  roll  and  pitch  system  result  from 
mass  unbalances  interacting  with  the  angular  scanning  motion  about  the  azimuth 
axis,  and  combine  to  form  a  single  spiral  torque  at  the  same  rate  and 
direction.  Its  effect  on  the  system  is  considered  in  Section  8.  Again  an 
infinite,  but  rapidly  converging  series  of  spiral  motions  are  produced,  as 
forcing  torques  are  matched  with  reaction  torques.  The  magnification  of  the 
forcing  torque  by  the  servo  loop  resonance  is  considered  in  terms  of  the 
torque  demand  from  the  roll  and  pitch  axis  motors. 

Section  9  considers  torques  in  the  azimuth  and  elevation  drives  under  the  same 
steady  conditions  postulated  for  the  preceding  sections. 

Section  10  considers  the  torques  about  all  four  motor-driven  axes,  due  to  a 
varying  elevation  angle  combined  with  otherwise  steady  scanning  conditions. 
Section  11  considers  the  same  torques  during  azimuth  axis  acceleration. 

If  the  azimuth  axis  is  not  truly  vertical,  due  to  roll  and  pitch  correction 
errors,  the  true  azimuth  and  elevation  angles  of  the  boresight  differ  slightly 
from  the  indicated  values.  The  differences  arc  evaluated  in  Section  12. 

Sections  13,  14,  15  and  16  consider  the  boresight  errors  and  changes  in  motor 
torque  demands  resulting  from  deviations  from  the  ideal  geometrical 
arrangements  of  the  four  motor  driven  axes. 

The  results  of  the  analysis,  presented  in  tabular  and  graphical  form,  provide 
a  firm  basis  for  the  design  of  the  servo  drives  and  of  the  procedure  and 
facilities  for  achieving  the  necessary  degree  of  mass  balance. 
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2.  IDEAL  ARRANGEMENT  OF  MOUNT  AXES 

In  this  Section  we  describe  and  define  the  ideal  arrangement  of  the  axes  and 
assemblies  of  the  elevation-over-azimuth-over-pitch-over-roll,  stabilised, 
four  axis  aerial  mount.  Deviations  from  the  ideal  arrangement  are  considered 
in  later  sections.  A  simple  extension  covers  the  case  of  two  aerials,  back  to 
back,  on  two  separate  elevatable  assemblies. 

The  basic  concept  of  the  ideal  arrangement  is  that  the  motor-driven  roll- 
correction,  pitch-correction  and  azimuth  axes  always  intersect  at  a  common 
point,  0,  irrespective  of  the  roll  and  pitch  correction  angles,  ti  and  ip . 
Furthermore,  the  pitch  and  roll  axes  are  always  perpendicular,  as  are  the 
pitch  and  azimuth  axes. 

The  roll  correction  axis  is  adjusted  to  be  horizontal,  and  parallel  to  the 
rear  to  front  centre  line  of  the  carrying  vehicle,  when  the  latter  is  in  its 
nominal  level  condition.  The  pitch  correction  axis  is  defined  to  be 
horizontal  under  the  same  vehicular  condition  and  when  the  roll  correction 
angle,  n ,  is  zero.  The  azimuth  axis  is  defined  to  be  vertical  (pointing 
towards  the  zenith)  when  the  vehicle  is  level,  ri  is  zero,  and  the  pitch 
correction  angle,  </> ,  is  zero. 

The  Cartesian  coordinate  system  with  the  b-suffix  (figure  1(a))  is  fixed  to 
the  base  of  the  aerial  mounting  system.  Its  origin  is  the  point,  0,  and  its 
x-axis  coincides  with  the  roll  correction  axis.  The  y  and  z-axes  coincide 
with  the  pitch  and  azimuth  axes  respectively  when  ti  and  >p  are  zero. 

The  Cartesian  system  with  the  g-suffix  is  fixed  to  the  gimbal  ring.  Its 
origin  is  also  0,  its  x-axis  also  coincides  with  the  roll  axis,  and  its  y-axis 
coincides  with  the  pitch  axis. 

The  Cartesian  system  with  the  s-suffix  is  fixed  to  the  stabilised  platform 
(figure  1(b)).  Its  origin  is  also  0,  its  y-axis  also  coincides  with  the  pitch 
axis,  and  its  z-axis  coincides  with  the  azimuth  axis. 

The  Cartesian  system  with  the  r-suffix  is  fixed  to  the  rotating  assembly.  Its 
origin  is  also  0,  its  z-axis  also  coincides  with  the  azimuth  axis 
(figure  1(c)),  and  its  y-axis  is  defined  to  be  parallel  to  the  motor  driven 
elevation  axis.  The  azimuth  angle,  0,  is  defined  to  be  zero  when  the  x  and  y- 
axes  coincide  with  those  of  the  s -system. 

The  elevation  axis  is  parallel  to  0Y^  (figure  1(d))  and  cuts  the  X^OZ^  plane 

at  the  point  0  ,  (Xl,  0,  Z ! )  in  r-system  coordinates.  The  Cartesian  system 

with  the  e-suffix  is  fixed  to  the  elevating  assembly.  Its  y-axis  coincides 
with  the  elevation  axis.  The  x-axis,  O^X^,  i-s  defined  to  be  the  boresight 

direction  of  the  associated  aerial.  The  elevation  angle,  0,  is  defined  to  be 
zero  when  the  bore  sight,  and  x-axis,  are  parallel  to  the  x-axis  in  the  r- 
system . 

The  Cartesian  system  fixed  to  the  second  (back  to  back)  elevating  assembly 
would  coincide  with  that  of  the  first  elevating  assembly  if  the  latter  were  to 
be  rotated  through  radians  about  the  0Z^  axis.  0^  would  occupy  the  point 

(-Xi,  0,  Zi)  in  the  r-system,  and  the  x  and  y-axes  would  be  oppositely 

directed  to  those  in  the  r-system.  An  increase  in  <f>  tilts  the  boresight  of 
the  second  assembly  upwards. 

We  define  the  following  rotation  matrices  and  position  column  vector. 
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3.  CALCULATION  AND  TRANSFORMATION  OF  TORQUE 
3.1  General  Cartesian  axes 


Consider  a  general  set  of  Cartesian  axes,  OXYZ,  moving  in  inertial  space, 
as  shown  in  figure  1(e).  Consider  also  a  point-mass,  of  magnitude,  m, 
located  at  the  time-dependent  coordinates,  x,  y,  z,  with  respect  to  these 
axes.  The  mass  exerts  a  reaction  force  as  a  result  of  its  absolute 
acceleration  and  of  its  presence  in  the  earth's  gravitational  field. 
Suppose  it  is  resolved  into  the  components,  F^,  F  ,  F^ ,  parallel  to  the  x, 

y  and  z-axes.  We  wish  to  determine  the  torque  components,  T^,  T^ ,  T^, 


about  the  same  axes. 


An  examination  of  the  situation,  using  figure  1(e),  shows  that  the 
following  matrix  relationship  applies. 
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L 

_ 

z_ 

We  define  the  column  vectors  of  position,  P,  force,  F,  and  torque,  T,  as 
well  as  the  matrix,  N(P),  where: 
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N 


P  =  [x  y  z  ] 


F  = 


[F  F  F  j 
x  y  z 


T  =  [ 


N(P)  = 


T  T  T  ] 
x  y  z 


[0  0  ll.P 


>  O) 


[0  0  -1] .P  [0  1  0] .P 

0  [-1  0  0] .P 


[0  -1  0] .P  [10  0  ]  .  P 


0 


Then  the  above  relationship  becomes: 


T  =  N(P).F 


(4) 


3.2  General  inertial  axes 

Consider  a  general  set  of  Cartesian  axes,  O'X'Y'Z',  fixed  in  inertial  space 
with  O'z'  pointing  towards  the  zenith.  Suppose  the  same  point-mass  is 
located  at  the  time-dependent  coordinates,  x1,  y '  ,  z' ,  with  respect  to 
these  axes.  Suppose  the  same  reaction  force  is  resolved  into  the 
components,  ' ,  F  ' ,  ' ,  parallel  to  the  corresponding  inertial  axes. 

We  define  the  column  vectors,  P'  and  F1 ,  and  L,  and  the  coordinate  rotation 
matrix,  R,  such  that: 


F'  =  [F  '  F  '  F  ’ 1T 
x  y  z 

P'  =  L  +  R.P 


Here,  L  defines  the  time-dependent  offset  between  the  origins  of  the  two 
Cartesian  systems,  and  R  defines  the  time-dependent  angles  between  the 
corresponding  pairs  of  axes. 

1 f  g  is  the  gravitational  acceleration,  then: 


F'  =  -  m. (G  +  P’) 

G  =[00  g]T 


(6) 


where  the  dots  denote  differentiation  with  respect  to  time. 

The  force-components  transform  isomorphical ly  with  the  position  coordinates 
under  rotation  of  axes,  and  are  unchanged  under  translation  of  axes. 
There  fore : 
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¥ '  =  R  .  F 


(7) 


Combining  the  relations  (4),  (6)  and  (7),  we  get  the  general  relationship: 


T  =  -m.N(P) .  R~  1 . (G  +  P* ) 


(8) 


3.3  Prope.rties  of  N(P) 

It  can  be  shown  (by  evaluation)  that  the  inverse  of  the  matrices,  A,  B  and 
C  (from  equation  (1))  are  obtained  by  changing  the  sign  of  the  angle 
parameter . 

It  can  also  be  shown  (by  evaluation)  that,  if  Q  is  a  matrix  of  the  type  A, 
B  or  C  ( ie  a  coordinate  rotation  matrix),  then  for  all  F: 


N(Q.P)  =  Q - N ( P ) .  Q  ~  1 
N  ( P ) .  P  =  0 


(9) 


3.4  Expansion  of  torque  relation 
I'sing  equations  (5)  and  (9): 

N(P).k'1.(Pi  -  I.)  =  R'1.N(R.P).idLt(Jt(R.P)} 

=  R'l-5-t{N(RP)-dt(RP)} 

=  ( R  ~  1  .  R  +  4r){N(P).(R~l.R.P  +  P)} 


Therefore  f rom  equation  (8): 

T  =  -m.(R*l.R  +  Ht){N(P).(R_I.R.P  +  P) }  -  m .  N(I>)  .  R_  1 .  (G  +  L)  (10) 


3.9  Axis  transformation 

We  wish  to  determine  the  torque  T" ,  given  T  and  the  axis  transformation 
re  1  a  t  i onsh  ip : 


parameters,  d  and  <;0- 


8.3  Magnification  of  motor  torque 

The  ratio  of  the  magnitude  of  the  servo  motor  torque,  to  that  of  the 
forcing  torque  (which  equals  the  total  deflection  induced  reaction  torque) 
is  of  interest  because  significant  magnification  can  occur  if  the  forcing 
angular  rate  approaches  the  servo  resonance  frequency.  We  consider  the 
simple  cast;  of  negligible  servo  unbalance  (u  and  f  are  zero)  and  one  spiral 
component  of  forcing  torque. 

The  expression  for  motor  torque  is  given  in  equation  (27).  If  values  for  U 
and  F  are  substituted  from  equation  (29),  and  for  (h+j-'/')  from 
equations  (31)  and  (42),  we  find  that  the  ratio  of  motor  torque  to  forcing 
torque  is  given  by: 

T  /T  =  -eSL.{w02  +  g.V/U+Q)  +  2c0  .w0  .s)/G(s)  (48) 


The  magnitude,  of  this  expression  is  plotted  in  figure  4,  parts  (a)  to  (d), 
for  fs  -  j.O),  as  a  function  of  w0/0,  and  for  several  values  of  the 

p  if  ime; ers ,  d  and  c0 •  V  is  assumed  to  be  zero.  For  (d  =  0.5)  the  ratio  is 

id..;  •  -v  a  !  1  y  vm  i  t.y  . 

hi*.'  ■  :  i’  : 


•  ■,  The  uiranete.r ,  V,  appears  in  equation  (48)  only  because  of  the 

it  ion  of  u0  in  equation  (29). 

For  values  of  d  greater  than  0.5,  the  system  is  in  stable 

•an; librium  even  without  the  motor  torque,  and  the  magnification  is  less 

t  i:an  an  i  ty  . 

'c.',  For  values  of  d  less  than  0.5,  the  system  is  unstable  without  the 
motor  torque,  and  magnifications  of  up  to  3  dB  can  occur. 

id)  if  the  servo  is  to  be  in  control  of  the  system,  and  particularly  if 

i  ; I  ,.ss  than  0.5,  the  value  of  w0  must  be  at  least  three  times  the 

v  ■!■'••  oi  0. 

9.  STEADY  STATE  TORQUES  IN  AZIMUTH  AND  ELEVATION 


[ n  i iiis  Section  we  consider  torques  about  the  azimuth  and  elevation  axes 
under  ' lie  same  steady-state  conditions  previously  assumed,  viz 
(0  0-0-0,  a  -  <1/  =  0).  The  assumptions  imply  that  these  torques  are 

resisted  by  the  respective  servo  motors  with  no  change  in  motion  (ie 
infinite  stiffness). 

r)  .  2  Azimuth  torque  (T  ) 

r  z 

We  use  the  combined  rotating  and  elevating  assembly  parameters  (Sub- 
sect  ion  4.5)  derived  from  fable  4,  to  f  ind  that: 
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(b)  For  the  H-term  from  equation  (40): 


s  ”  j  .  0 


Tf  =  j.H.g.e- 


(c)  For  the  K-term  from  equation  (35),  excited  by  the  spiral  deflection 
defined  in  equation  (31): 


s'  =  s  +  2j  .  8 


Tf  '  =  -  K .  A  .  e”-^  .  {  (s  )  2  +  2  j  .  0  .  s  "} 


(d)  For  the  K-term  from  equation  (35),  excited  by  the  main  spiral 
deflection  produced  by  the  forcing  torque  T^.: 


s'  =  s"  +  2j.G 


(Tf ' )" 


K.Tf .eje. (s2-2j.0.s)/{(J+Q).G(s)} 


(A')"  =  K.Tf  .ejP.  (2j  ,0.s-s2)/{  (J+Q)2  ,G(s)  .G:(s’"’+2j  .0)} 


Note  that,  if  s  =  j.O,  then: 


s  -  j.O 


Tf  -  j  .  L .  0 2  .  e~ 


s'  =  j.O 


(2j.0.s-s2)  =  -  0‘ 


(s  +2 j . 0 )  =  j.O 


and  the  output  angular  frequency  and  direction  are  the  same  as  the 
input.  This  does  not  apply  for  (s  =  -j.O),  for  which  the  output  occurs 
at  (s ' =  3 j . 0 )  . 

Figure  3,  parts,  (a)  to  (e)  ,  show  the  value  of  the  magnitude  of  G(s) 
(normalised  to  02)  plotted  as  a  function  of  w0  (normalised  to  0)  for 


=  - 1.01  a 


iMQ 


i.0).  and  for  several  values  of  the 
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7.3  Perturbations 

Newton's  approximation,  or  some  other  iterative  method,  can  be  used  to 
correct  the  value  of  s,  determined  in  Sub-section  7.2,  to  make  the 
function  (33)  zero  for  finite  values  of  u  and  f.  This  converges  rapidly 
for  values  of  the  order  of  0.2  or  less. 

If  the  value  of  K  is  small  but.  finite,  a  spiral  torque  can  be  determined 
from  equation  (35)  for  each  of  the  four  spiral  deflection  terms  (two  values 
of  s,  and  their  conjugates).  These  are  then  regarded  as  forming  torques, 
which  produce  spiral  deflection  pairs  at  the  same  angular  frequency,  as 
discussed  in  Section  8.  These  in  turn  produce  further  reactions  from  the 
linear  reaction  section.  Normally  the  infinite  sequence  of  terms  converges 
rapidly  to  negligible  amplitude.  (Note  that  the  terms  from  equation  (35) 
need  a  sign  reversal,  to  convert  them  to  forcing  torques). 


8.  STABILISATION  SYSTEM  FORCED  RESPONSE 

8 . 1  General 

An  examination  of  Table  6  shows  that,  under  the  same  conditions  assumed  in 
Sections  6  and  7  (ie  0  =  0  =  0  =  0,  and  ri  -  ip  -  0),  there  are  still  some 
torques  not  accounted  for  as  deflection  induced  torques  (Section  6)  .  These 
forcing  torques  are  independent  of  any  (small)  spiral  deflection  motions, 
and  are  given  by: 

(T  ,,+j.T  .  ,)  =  j.L.02.ej9l  +  j.H.g.ej0h  (40) 

roll  J  pitch  J  jo  v 

Those  can  be  generalised  to  the  form: 


(Troll+j •Tpitch) 


where,  in  this  case,  (s  =  j.O)  and  T^  is  the  appropriate  complex  quantity 
for  equivalence  of  equations  (40)  and  (41). 

The  forcing  torque  must  always  be  balanced  by  a  motion  induced  reaction 
torque  (from  equation  (33)).  Thus  s  is  defined  by  the  forcing  torque,  and: 


(J+Q).{G(s)  -  f(u,f)} 


(42) 


where  f(u,f)  represents  the  right-most  term  in  equation  (33). 

The  induced  spiral  deflection  terms  are  then  defined  by  equation  (32).  The 
main  term  is  clearly  at  the  same  angular  rate,  and  in  the  same  direction  as 
the  forcing  torque. 

8.2  Special  cases 

(a)  For  the  L-torm  from  equation  (40): 
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This  is  applicable  independent,  ly  to  the  roll  and  pitch  loops, 
because  there  is  no  cross  coupling.  Critical  damping  occurs  with 
(cr  1). 

(b)  If  0  f  0,  we  substitute  (s  =  a  +  j.w)  in  (36)  and  equate  the  real 
and  imaginary  parts  separately  to  zero.  Thus: 

(w-d. 0) . (a+c0 . w0 )  +  c0.w0.d.8  =  0 

-(w-d.O)2  +  (a+c0.w0)2  +  d2.02  +  w02 . (l-c02)  =  0 

Note  that  (w  =  d.0)  is  not  a  solution,  else  (c0.w0  =  0)  and  o  is 
imaginary.  Therefore: 

a  =  -c0  - w0  •  <*>/ (wd •  9)  (38) 

and,  multiplying  by  -(w-d.8)2  and  rearranging: 

2 (w-d . 0 ) 2  =  d2 . 0  2  +  w02. (l-c02) 

+  {  (d2  .  0  2  +  u>o2.(l-c02))2  +  4c0  2 . w0  2 ■ d2 . 02 
=  2d2 .0Z  +  w02  . (l-c02)  -  d2  .  02 

+  {(w,2.(l-co2)-d2.02)2  +  4wo2.d2.02}*  (39) 

Note  that,  from  equations  (39), (38): 

(1)  (w-d . 0 ) 2  >  d2 .G2 

Thus:  w  >  2d.0  or  w  <  0 

(2)  In  both  cases,  a  <  0,  and  transients  decay  to  zero  with  increase 
in  time. 

(3)  I  — I  =  --~d' ?  I  takes  the  same  value  for  both  solutions,  so  the 

I  a  |  c  0.oj0  I 

spirals  are  symmetrical,  although  one  is  executed  much  faster  than 
the  other. 

Figure  2  shows  the  magnitude  of  (w  -  d.0). from  equation  (39),  plotted  in 
a  normalised  form,  as  a  function  of  wo/d.0  and  c0.  It  can  be  seen  that, 

for  large  values  of  c0 ,  there  is  a  forward  spiral  response  at  an  angular 

rate  of  just  greater  than  2  d.0,  together  with  a  reverse  spiral  resoonse 
at  an  angular  rate  close  to  zero.  The  directions  are  unchanged  but  the 
rates  increase  considerably  as  the  damping  factor,  c#  ,  is  reduced, 

depending  also  on  the  stiffness  of  the  servo  loops  (w0)- 
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Then  equation  (30)  becomes: 


'(Troirj-Tpitch)/(J+Q) 


st 


A.e  ,{G(s)  -  (u.w02+  2f.c0.u0-s)V(G(s)+  4j.d.0.s)}  (33) 


6.4  Non-linear  reaction  to  spiral  deflection 

We  now  consider  the  contribution  of  the  K-term  to 
Table  6: 


the  spiral 


torque . 


From 


(Troll+j ■Tpitch) 


K.  e2j  ‘  9k.  {  C'ri-j  .'^)  +  2j  .  0  .  (n- j  .^) } 


=  K.A*.e(s  t+2j  ■9k)  .{(s'")2  +  2j.0.s*} 


under  the  single  spiral  deflection  (equation  (31)). 

The  significance  of  this  component  will  be  considered  later. 


7.  STABILISATION  SYSTEM  TRANSIENT  RESPONSE 


(34) 


(35) 


7 .  1  Genera  1 

The  transient  response  of  any  system  is  obtained  by  equating  to  zero,  the 
total  deflection-induced  torque.  In  this  case  it  is  the  sum  of  the 
functions  (30)  and  (34).  A  closed  solution  can  not  be  obtained  directly, 
However  we  can  obtain  a  solution  by  assuming  that  K,  u  and  f  are  zero,  and 
then  perturbing  it  to  take  into  account  the  small  finite  values  of  these 
parameters . 

7.2  Unperturbed  solution 

We  assume  the  form  of  the  solution  is  as  defined  in  equation  (32).  Then, 
with  (K  =  u  =  f  =  0) ,  from  equations  (33)  and  (35),  the  torque  is  zero  if: 


G(s)  =  0 


This  can  be  rearranged  algebraically  to  give: 


(s- j .d. 0+co .w„ )2  +  2 j . c0 . w0 . d . 0  +  d2.02  +  u02.(l-c02)  =  0 


(36) 


(a)  if  0  =0,  wc  get  the  well  known  result: 


=  -c0  .  w0  ±  u0 


( c  o  ‘ 


-1) 


i 


S 


(37) 
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U-8-V 
J  +  Q 


2{(U-g.V).(J+Q)}‘ 


d 


and  assuming  that  — -  0,  equation  (28)  becomes: 


'(Troll+j  •Tpitch)/(J+Q)  =  +  2(co.w,-j.d.0).(n+j.*) 


+  Wo2  ■  (h+j  -<f>) 


+  u.u02 .  (n-j  .\i>)  +  2f . c0 . w0 .  (n •  j  .<f>) 


6.3  Linear  reaction  to  spiral  deflection 


Suppose : 


n+j.if>  =  A.e~ 


s  =  a+j.w 


From  parts  (a)  and  (b)  of  figure  1,  this  is  an  anti-clockwise  spiral  motion 
of  the  point  Z (looking  down  towards  the  origin)  about  the  point,  Z^,  with 

angular  frequency,  u>,  and  logarithmic  amplitude  increment,  o.  A  should  be 
considered  as  complex,  to  represent  the  amplitude  and  phase  angle  at  t  =  0. 
Replacing  s  with  its  complex  conjugate  reverses  the  direction  of  the 
spiral . 

If  we  substitute  equation  (31)  into  equation  (30)  we  obtain  a  major  term  in 

eSt  (being  a  spiral  torque  in  the  same  direction  as  the  spiral  deflection, 

* 

,  s  t 

and  with  the  same  angular  frequency  and  increment)  and  a  minor  term  in  e 
(being  a  spiral  torque  in  the  reverse  direction  with  the  same  angular 
frequency  and  increment)  due  to  the  unbalanced  servo  motor  gains  u  and  f. 
It  is  more  useful  to  add  a  small  spiral  deflection  in  the  reverse  direction 
to  cancel  the  reverse  spiral  torque.  This  is  achieved  if: 


(h+j-iJO  = 


,  St  *  s  t 
A  .  e  -  B  .  e 


G(s)  =  s2  +  2 (c0 . w0 - j . d . 9) . s  +  we 


B  =  A.(u.w02  +  2f . c0 . u0 • s)/ (G(s )  +  4j.d.8.s) 
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2K 

e 


m.(z  V 
o 


m.  (x  '  ): 
e 


(26) 


Thus  K  is  proportional  to  the  difference  between  the  principal  moments  of 

inertia  in  the  X  Z  -  plane.  For  =  0,  the  two  principal  moments  of 

inertia  must  be  equal  (and  the  e' -system  axes  are  principal  for  all  values 
of  <p )  . 


6.  DEFLECTION  INDUCED  TORQUES  IN  THE  STABILISATION  SYSTEM 
6 .  1  Genera  1 

Tn  this  Section  ue  examine  the  torques  present  in  the  roll  and  pitch 
correction  system  due  to  the  deflection  of  the  azimuth  axis  from  the 
vertical  position  (n+j-<i>)  and  its  time  derivatives,  under  the  assumptions 
(0  =  0  =  0  =  0)  and  (n  =  ip  -  0).  These  consist  of  gravitational  and 
inertial  reactions  (already  determined,  and  listed  in  Table  6)  and  servo 
motor  supplied  torques. 

We  assume  that  the  roll  and  pitch  correction  servo  motors  produce  torques, 
proportional  to  the  deflection  error  and  its  first  time  derivative,  and  in 
the  direction  to  reduce  that  error.  We  wish  also  to  examine  the  effect  of 
small  differences  in  the  constants  of  proportion,  between  the  roll  and 
pitch  systems.  We  therefore  postulate  the  following  relations  for  motor 
torque : 


T 

roll 

=  -  U.  (1+U).u  -  F.  (l+f).r| 

\ 

T  . 
p  l  tch 

=  -  U.  (1-u)  .ip  -  F.  (1-f )  .tp 

> 

(27) 

fTroll  +  j 'Tpitch^ 

=  u.  (n+j  .<i>)  +  F.  (n+j  . 0) 

+  u.u.  (n-j  -<p)  +  F.  f.  (n-j  -V>) 

/ 

6.2  Linear  reaction  torque 

From  Table  6,  and  the  assumption  (0  =  0)  all  reaction  terms  are  linear 
except  those  in  K.  For  this  sub-section  we  assume  that  K  =  0 .  Then  the 
total  reaction  torque  is  given  by: 


-IT  ,,+  j.T  .  .  )  -  (  HQ)  -  (n  +  j  •'</-)  +  (F-2  j  ..1.0  )  .  (n+j  .<H 

roll  pitch 

+  d'-g.  V) .  (n+j  .</>) 

*  ii  c  .  fn-  j  .<H  +  f .  F .  ( n  -  j  )  (28) 


Making  the  normalising  sebsi  i  '  ■  i *  u>ns: 
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(2)  The  conditions  involving  H  and  H  are  the  obvious  ones  that  the 

centre  of  mass  of  the  elevating  assemblies  be  on  the  elevation  axis,  and 
that  of  the  rotating  plus  both  elevating  assemblies  be  on  the  azimuth 
axis.  The  structure  is  then  statically  balanced.  These  conditions  can 
be  easily  met  in  practice,  with  proper  mechanical  design  and  with 
adjustable  weights  for  fine  trimming. 

(3)  The  conditions  involving  and  are  those  for  dynamic  balance  of 

the  elevating,  and  total  rotating  assemblies,  about  the  elevation  and 
azimuth  axes,  respectively.  These  conditions  are  less  easily  met,  and 
require  careful  attention. 

(4)  The  condition  involving  is  very  difficult  to  achieve  in 

practice.  Its  phy..ical  implication  is  discussed  in  detail  in  Sub¬ 
section  5.6.  With  two  elevating  assemblies,  the  effect  about  the  roll 
and  pitch  axes  tends  to  cancel,  if  the  two  values  of  <f>  are  not  too 
dissimilar.  This  may  not  be  a  reasonable  constraint,  however. 

(5)  If  torque  balance  is  required  for  small  constant  values  of  ip  and  ri 
(as  well  as  at  the  zero  values)  then  the  Tables  6,  1,  3,  4  and  2  show 
the  extra  requirement  to  be: 


V  +  M  . Z !  =  0  (25) 

r  e 

The  e-suffix  term  implies  the  difference  between  the  two  elevating 
assemblies . 

5.6  Physical  significance  of 

The  Cartesian  system  with  the  e' -suffix  is  defined  by: 


P  '  =  D.P 

e  e 


Then,  by  evaluation: 


m  •  (x  ' ) 2 

G 


C 


L 


m.(z 

e 


/ 

e 


ID.X  '  .  Z  * 

e  e 


J  -  K  . cos (2$  -  3  ) 
e  e  e 


J  +  K  .cos(20  -  3  ) 
e  e  e 


-  K  .sin(2 <t>  -  3  ) 
e  e 


We  choose  20  =  3  ,  making  the  cross  product  function,  zero.  Hence  the 
e 

e* -system  becomes  a  principal  set  of  axes,  and  the  other  two  functions  take 
the  maximum  and  minimum  possible  values.  Then: 
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5.2  Roll  and  pitch  torques 


(T  +  j.T  .  .  )  =  j.02.ej°.{(L  .cos*  -  K  .sin20, 

roll  pitch  r  r  e  k 


+  M  .Xj.Zi  +  Xj.H  .  sin0,  +  Zx-H  .  cos0, ) 
e  e  h  e  h 


+  j.(L  .sin?  +  L  .sin  0,  +  Zx.V  )} 
r  r  e  1  1  e 


JQ 


+  j.g.eJ  .{(H^.cosa^  +  He>cos0^  +  M^Xx) 


+  j  .  (H  .  s  ina  +  V  )  } 
r  r  e 


(21) 


For  this  expression  to  be  identically  zero  for  arbitrary  values  of  0,  the 
requirements  are: 


H  =  K  =  L  =  0 

e  e  e 

H  .cosa  +  M  .Xi  =  0 

r  re 

H  .sina  +  V  =  0 

r  re 

L  .cosT  +  M  .Xi.Zj  =  0 

r  r  e 

L  .siny  +  Zx.V  =  0 

r  r  e 


\ 


>  (22) 


5.3  Azimuth  torque  (T  ) 

rz 

This  torque  is  automatically  zero,  so  no  special  requirements  apply. 

5.4  Elevation  torque  (T  ) 

T  ,  =  GZ.{X1.H  .sin0,  -  K  ,sin20  }  +  g.H  .cos0,  (23) 

elev  e  he  k  e  h 

For  this  expression  to  be  identically  zero  for  arbitrary  values  of  0,  the 
requirements  are: 

H  =  K  =  0  (24) 

e  e 


5.5  Discussion  of  results 

(1)  The  conditions  for  e-suffix  terms  in  equation  (22)  imply  the 

1 9 

difference  between  the  two  elevating  assemblies,  because  eJ  terms  tend 

to  cancel  (Sub-section  4.6).  But  if  different  values  of  0  are  allowed 

for  the  two  units,  cancellation  is  not  possible  and  the  zero  condition 

for  H  ,  K  and  L  apply  to  each  unit  individually.  The  conditions  in 
e  e  e 

equation  (24)  apply  individually  in  either  case. 
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T 


roll 


T  =  [100  I.B.C.T 

gx  r 


T 

P 


itch 


T 

sy 


[010] .C.T 


Both  T  ,  .  and  T  .  ,  are  real  variables  because  B,C  and  T  are  real.  (Note 

roll  pitch  r 

that  complex  terms  in  occur  as  the  sum  of  conjugate  pairs).  Therefore: 


\ 


T 

roll 

— 

Re{Troll+  j- 

T  .  ,  } 

pitch 

T 

pitch 

= 

Im{Troll+  j- 

T  .  u} 
pitch 

>  (20) 

(T 

roll  +  ^ ‘ Tpitch^ 

= 

[e-^0  j  .e-^0  ip 

].T 

r 

/ 

Table  5  lists 

the  contribution 

of  the 

elevating 

assembly  to 

(Troll+j •Tpitch) ’ 

evaluated  from 

Table  3,  using 

equation 

(20). 

Table  6 

lists  the  contribution  of  the  combined  assemblies,  at  fixed  elevation,  and 
is  derived  from  Table  4  by  deleting  the  r-suffix  from  the  mass  distribution 
parameters  (as  discussed  in  Sub-section  4.5). 

5.  CONDITIONS  FOR  BALANCE  WITH  STEADY  OPERATION 
5.1  Steady  operation 

For  a  continuously  scanning  radar,  it  is  normally  very  desirable  that  all 
torques  about  the  four  motor  driven  axes  (roll,  pitch,  azimuth  and 
elevation)  should  be  zero  when: 

(a)  the  aerial  mount  is  fixed  to  a  level  surface 

(b)  the  rotating  assembly  is  fully  erect  (ie  if>  =  n  =  0) 

(c)  there  is  a  fixed,  but  arbitrary  elevation  angle  (0) 

(d)  there  is  a  constant  (or  zero)  scan  rate  (9). 

Under  these  conditions,  the  torques,  obtained  from  Tables  6,  1,  3,  4  and  2 
become  as  shown  in  the  next  sub-sections. 
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of  torque  in  the  9-column  requires  that: 

L.coslf  =  L  .cos2f  +  M  .Xj.Z!  -  K  .sin20,  +  H  . (Xi.sintf,  +  Z^cos^,  ) 
r  r  6  6  K  6  it  ri 

All  such  relationships  are  listed  in  Table  1. 

Note  that  most  of  the  non-suf fixed  parameters  are  functions  of  <fi . 

4.6  Contribution  of  second  elevating  assembly 

As  discussed  in  Section  2,  the  second  elevating  assembly  is  (ideally) 
equivalent  to  the  first  assembly  rotated  by  n  radians  about  the  z-axis  of 
the  r-system.  Its  contribution  to  its  own  is  therefore  given  by  the 

terms  in  Table  2,  except  that  each  0  (including  any  in  the  column  heading) 
is  replaced  by  (0+ff) .  The  effect  of  this  is  to  change  the  sign  of  entries 
in  those  columns  containing  exp(+j.0)  in  the  heading. 

The  contribution  to  T^  is  determined  indirectly.  Clearly,  the  contribution 
to  Tg  is  obtained  by  replacing  0  by  (0+Jr)  wherever  it  occurs.  Thus  the 
contribution  to  T^  is  given  by: 


(19) 


where  the  prime  denotes  replacement  of  0  by  (0+7T).  Thus  the  rule  for 
Table  3  becomes:  change  the  sign  of  elements  in  the  x  and  y  components  of 
torque;  then  change  the  sign  of  elements  in  the  columns  containing 
exp(+j.0)  in  the  heading. 

Note  that  this  algorithm  for  T^  is  equivalent  to  changing  the  sign  for  each 

occurrence  (in  element  and  row  and  column  heading)  of  H  ,  V  ,  Xi  ,  $,  , 

e  e  n 

0^.  Executing  this  latter  algorithm  on  the  relationships  for  the  non- 

suf  fixed  mass  distribution  parameters  in  Table  1,  results  in  a  change  of 
sign  for  contributions  to  those  parameters  which  change  sign  when  both  Xf 

and  y^  are  changed  in  sign,  ie  H.cosa,  H.sina,  L.cos2T  and  L.sinJf.  This  is 

as  would  be  expected,  since  the  x  and  y  directions  in  the  r-system  are 
reversed  for  the  second  elevating  assembly. 

Final  sign  reversals  for  contributions  to  T  and  Tg  reflect  the  balancing 

effect  of  the  second  assembly  upon  the  first,  due  to  the  symmetrical 
structure  and  location  about  the  azimuth  axis. 

4.7  Torque  about  roll  and  pitch  axes 

From  figure  1  and  equation  (12)  it  follows  that: 
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From  equations  (2)  and  (11): 

K  =  E 
Q  =  D 


Items  (k)  to  (u)  of  Appendix  I  are  additional  key  results  in  the  evaluation 
of  equation  (16).  The  final  results  are  shown  in  Table  3  (in  the  same 
format  as  used  in  Table  2),  except  that  terms  in  (ie  those  which 

consider  the  elevating  assembly  as  a  point  mass  at  its  origin)  have  been 
included  in  Table  4,  together  with  the  contributions  of  the  rotating 
assembly . 

Note  that  the  z-component  of  T  is  the  torque  (due  to  the  elevating 
assembly)  to  be  resisted  by  the  azimuth  servo  motor. 

4.4  Contribution  of  rotating  assembly  to  T^ 

The  torques  generated  about  the  axes  of  the  r-system  due  to  the  reactions 
from  the  masses  forming  the  rotating  (but  not  elevating)  assembly  are  given 
by  equation  (14),  where  in  this  case  from  equations  (2)  and  (5): 

P '  =  P.  =  A.B.C.P 
b  r 


L  =  0 

R  =  A.B.C 


Items  (v)  to  (y)  of  Appendix  I  are  additional  key  results  in  the 

evaluation,  which  also  uses  the  set  of  ten  mass  distribution  parameters  for 

the  rotating  assembly  (indicated  by  the  r-suffix)  listed  in  Table  1.  The 

final  results  are  shown  in  Table  4,  which  includes  in  addition,  the  terms 

in  M  from  the  elevating  assembly,  derived  in  Sub-section  4.3. 
e 

Note  that  columns  involving  <f>  and  its  derivatives  have  been  deleted  from 
the  table.  The  mass  distribution  of  the  rotating  assembly  (alone)  is 
obviously  independent  of  the  elevation  angle  of  the  elevating  assembly. 

Note  also  that  the  z-component  of  T  is  the  torque  (due  to  the  rotating 
assembly  alone)  to  be  resisted  by  the  azimuth  servo  motor. 

4.5  Combined  contribution  to  T 

r 

It  is  obvious  that,  at  a  fixed  elevation  angle  (0),  the  elevating  assembly 
can  be  considered  as  an  extension  of  the  rotating  assembly.  Therefore  we 
can  define  a  combined  set  of  mass  distribution  parameters,  M,  J,  Q,  K,  L, 
H,  V,  a,  B  and  %,  which  behave  isomorphical ly  with  the  r-suffixed 
equivalents,  but  which  include  the  contribution  from  the  elevating  assembly 
at  a  fixed,  but  arbitrary,  elevation. 

The  relationship  between  the  non-suffixed  parameters  and  the  r  and  e- 
suf fixed  ones  can  be  determined  from  the  data  in  Tables  3  and  4,  by 
equating  total  torques  in  each  axis-column.  For  example,  the  x-component 
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where  from  equations  (2)  and  (5): 

P'  =  Pb  =  A. B.C. (E  +  D.P  ) 

P  =  P 

e 

L  =  A.B.C.E 

R  =  A.B.C.D 


(15) 


This  can  be  evaluated  directly,  using  the  mass-distribution  parameters 
defined  in  Table  1.  Between  them,  the  ten  parameters  associated  with  the 
elevating  assembly  (indicated  by  the  e-suffix)  cover  the  entire  range  of 
zero  order,  first  order  and  second  order  moments  of  mass,  measured  in  the 
coordinate  system  fixed  to  that  assembly.  The  form  of  the  parameters  has 
been  chosen  to  suit  the  evaluation  of  the  torques. 

Items  (a)  to  (j)  of  Appendix  I  are  key  results  in  the  step  by  step 
evaluation  of  equation  (14).  The  final  results  are  presented,  in  Table  2, 
in  the  form  of  a  matrix.  Each  entry  in  the  matrix  is  to  be  multiplied  by 
its  row  and  column  headers,  before  being  accumulated  into  the  total 
component  of  the  torque. 

The  first  five  columns  contain  dynamic  entries  due  to  the  motions  of 
rotation  (0)  and  elevation  (<t>) .  The  ninth  and  tenth  columns  contain  static 
entries  due  to  gravity.  Columns  6,  7  and  8  contain  dynamic  entries  due  to 
the  changing  values  of  roll  (n)  and  pitch  (tp)  correction.  Note  that  each 
of  columns  6,  7,  8  and  10,  represent  two  separate  sets  of  entries,  denoted 
by  the  upper  and  lower  sign  (where  a  choice  is  provided).  The  upper  sign 
is  used  for  one  set;  the  lower  sign  for  the  other.  Wherever  only  one  sign 
is  given,  it  applies  to  both  sets. 

The  complex  variable  form  of  column  headings  simplifies  evaluation  and 
presentation,  and  will  be  found  to  be  useful  in  later  Sections. 

Note  that  the  y-component  of  Tg  (as  listed  in  Table  2)  is  the  torque  to  be 
resisted  by  the  elevation  servo  motor. 

4.3  Contribution  of  elevating  assembly  to  T^ 

Consider  the  torques  generated  by  the  reactions  of  the  masses  located  on 
the  elevating  assembly,  about  the  axes  of  the  r-system.  Evaluation  is 
easier  if  we  use  the  transformation  relationship  (12),  rather  than  the 
direct  relationship  (10).  Thus: 

T  =  Q.T  -  Y  m.N(K).Q.R'1.(G  +  P')  (16) 

r  e  L 

e 


where  the  variables  are  defined  in  equations  (15)  and  (17). 
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P'  =  L  +  R .q" 1 . (P"  -  K) 

and  from  equation  (8): 


T"  =  -  m.N(P") . (R.Q'1)"1 . (G  +  P’) 

=  -  m.N(K  +  Q. P) .Q.R* 1 . (G  +  P') 


=  Q.T  -  m.N(K) .Q.R- 1 . (G  +  P' ) 


(12) 


Note  that  the  second  term  of  equation  (12)  contains  no  second  order  terms 
in  the  elements  of  a  position  vector.  All  second  order  terms  are 
transformed  within  the  first  term  of  equation  (12). 


4.  TORQUES  DUE  TO  MASS  UNBALANCE 


4.1  General 


In  this  Section  we  consider  the  torques  acting  about  various  axes,  with  the 
aerial  base  levelled  on  firm  ground,  and  with  the  azimuth  axis 
approximately  vertical.  In  this  case,  the  roll  and  pitch  correction 
angles,  ti  and  4>,  are  small  enough  that  second  order  terms  in  these 
quantities  can  be  neglected.  Then: 

A .  B  -  I  +  X  I 


X 


r  0  0  i p- 

o  o  -n 
--<l>  n  o  _ 


(A  .  B )  '  1  =  I  -  X 


S  (13) 


dt 


(A .  B) 


X 


The  torques  considered  are  caused  by  mass  unbalances  within  the  structure 
of  the  rotating  and  elevating  assemblies.  To  calculate  them,  we  must  sum 
the  effect  of  reaction  forces  over  the  entire  rotating  and  elevating 
assemb  1  ies . 

On  the  assumption  that  the  b-system  of  axes  is  inertial,  with  the  z-axis 
pointing  towards  the  zenith,  this  system  is  isomorphic  with  the  singly 
primed  inertial  system  discussed  in  Section  3,  and  the  torques  can  be 
evaluated  using  the  relationships  derived  therein. 

4.2  Contribution  of  elevating  assembly  to  T^ 

Consider  the  torques  generated  about  the  axes  of  the  e-system  due  to 
reactions  from  the  masses  forming  the  normal  elevating  assembly.  From 
equation  (10): 
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L.{e  '  0  1  .  (n+j  .  V)  +  .  Cn-j  vt>")) 


-  £g.H.{e  ^  .  (T|+j  .  #)  +  e^0h.  ( n  -  j -<P)  } 


Re{L.e  ^  9  1  .  (ri+j  .#)  -  g.H.e  ^  9h  .  (n+j  •  4>) } 


where  Re  denotes  "the  real  part  of". 

Considering  a  single  spiral  deflection  component,  as  defined  in 
equation  (31): 

T  =  |  A  |  .Re{L.s2.e*St‘j(er“A)}  -  g  .  H  .  e{st‘ j  (0h"AA)  }  }  (50) 


In  the  common  case  (s  =  j.0): 


T  =  -  |  A  |  . L . 0  2 . cos  (AA  -X)  -  |  A  |  .  g .  H  .  cos  (AA-a)  (51) 

Both  terms  consist  of  a  steady  torque,  of  magnitude  dependent  on  the  phase 
relationship  between  the  relevant  mass -unbalance  and  the  deflection 
component . 

If  (s  =  -j.0)  or  (s  =  3 j . 0 )  the  torque  components  are  sinusoidal,  and  at  an 
angular  frequency  of  20. 

9.3  Elevation  torque  (T  ) 

From  Table  2,  and  assuming  the  same  single  spiral  deflection 

(equation  (31)  )  : 

T  ,  =  02.(X,.H  .sin#,-K  .sin2#,  )  +  g.H  .cos#, 

elev  1  e  he  k'°e  h 

+  Re[A.e^St  ^  0  ^ .  {  j  .  s 2  .  (2J  +Xj.H  .cos#  +  Zj.H  .sin#,  -j-L  .cos#.) 

6  S  il  6  u  6  I 

+  20.s.(J  +Xt.H  .cos#  -K  .cos2#  -j.L  .cos#  )-j.g.H  .sin#  }]  (52) 

6606  K  6  1  6  n 


For  (s  =  j . 9) : 


T  .  =  02.(Xj.H  .sin#,  -K  .sin2#,  )  +  g.H  .cos#, 

elev  e  he  k  e  h 


+  02. |A| . (Zj.H  .sin#,-X1.H  .cos#  +2K  .cos2#  ).sin^A 
g  n  g  n  g  1C 


+  0 2 . | A | . Lp . cos# ^ . cosAA  +  | A | . g . . s in#u . sinAA 


e  *V 


For  (s  =  j.w),  the  term  in  becomes: 


(53) 
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2  |  A  |  .  w .  (w-Q).sin(w.  t-Q+LA) .  J  (.54 ) 

Note  that,  for  the  second  elevating  assembly,  0  is  incremented  by  it  in 
equation  (52),  and  all  terms  in  A  take  a  reversed  sign. 

10.  TORQUES  DUE  TO  VARYING  ELEVATION  ANGLE 

Tables  5,  3  and  2  list  all  significant  torques  due  to  #  and  #,  about  the  roll 
and  pitch,  azimuth  and  elevation  axes,  respectively.  Major  terms,  and  those 
of  first  order  of  smallness,  are  given  below. 

(T  . .+  j.T  .. 
roll  pitch) 

i  fi  •  • 

=  j.eJ  .{<#».( 2 J  +Xi-H  .cos#,+Zj.H  .sin#,+j.L  .cos#  ) 
eenenei 

-  #1.(X1.H  .sin#,  -Zj.H  ,cos#,+j.L  .sin#.) 

e  n  e  n  e  1 

+  2j.0.#.(J  +Zi.H  . sin#. +K  .cos2#  ))  (55) 

6  6  il  6  K 


T  =  #.L  .sin#,-#2.L  .cos#, 

az  e  1  e  i 


+  20 . # . (X j . H  .sin#  -K  ,sin2#  ) 

6  fl  6  K 


+  2#.  |A|  .Je.Re[s.e(st"j0+jZjA)J 


T  .  =  2J  .# 

elev  e 


Note  that  all  terms  in  equation  (55)  tend  to  cancel  when  two  elevating 
assemblies  have  the  same  variation  of  elevation  angle.  Only  the  term  in  |A| 
does  so,  in  equation  (56). 


11.  TORQUES  DUE  TO  AZIMUTH  ACCELERATION 

Tables  6,  4  and  2  list  the  significant  torques  due  to  0.  Major  terms,  and 
those  of  first  order  of  smallness,  are  given  below. 
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T  =  -  2J.0 

az 


(59) 


8.L  .sin$, 
e  1 


12.  EFFECT  OF  ROLL  AND  PITCH  ERRORS  ON  BORESIGHT 
12.1  Bores ight  error 

As  discussed  in  Section  2,  the  aerial  boresight  is  assumed  to  be  the  axis, 
0^  X  ,  which  is  normal  to  the  elevation  axis  at  point,  0g,  and  parallel  to 

the  XY-plane  of  the  r-system  when  <f>  =  0  (see  figure  1(d)). 

Suppose  the  boresight  passes  through  a  fixed  target,  located  at  the  point, 
P  ,  under  two  separate  conditions.  Firstly,  with  no  roll  and  pitch  errors: 

A.B  =  I 

C  =  C 

D  =  D 


Secondly,  with  finite  roll  and  pitch  errors: 

A.B  =  A.B 

C  =  C  +  AC 

D  =  D  +  AD 

P  =  (r  +  Ar) .Z 
e 

Using  equation  (2)  to  determine  P^  for  each  case,  and  then  equating  che 
results : 


A.B. (C+AC) . (E  +  (D+AD) . (r+Ar) .Z) 


C . (E  +  D.r.Z) 
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Therefore : 


D_  1 .  ( I  +  C ~ 1 . AC ) . D . { ( I  +  D'1 .AD) . (r+Ar) .Z  +  D_1.E> 

=  D_ 1 . C~ 1 . (A . B ) ” 1 . C . D . ( r . Z  +  D_1.E)  (61) 

Assuming  that  Ar/r  and  the  elements  of  D  '.E/r  are  negligibly  small 
compared  with  unity,  then  equation  (61)  reduces  to: 


D‘ 1 .C*1 .AC.D.Z  +  D_ 1 . AD . Z  +  D_ 1 .  C~ 1 . AC . D . D~ 1 . AD . Z 
=  D"1.C*l.{(A.B)"1-I}.C.D.Z 


(62) 


Noting  that: 


cos (0+A0) 

=  cosA0.cos0  +  s inA0 . ^ (cos0) 

s in (0+A0 ) 

=  cosAO.sinO  +  si  iA0  .^(sin0) 

C  +  AC 

=  cosAO.C  +  sinAO.j^ 

do 

C~ 1 . AC 

j  p 

=  (cosA8-l).I  +  sinAO.C  l.— 

do 

and  using  equation  (1),  the  left  hand  side  of  equation  (62)  becomes: 


cosA8.cosA0  -  1 
s inA0 . cos (0+A0) 
sinA0 . cosA8 


(63) 


Making  the  substitutions: 


T1  =  JJ.COS0  1 


tp  =  y.sinfl 


>  (64) 


v  =  (n2+  1 1>2)^ 


and  including  terms  of  second  order  in  q  and  <p  from  relation  (1),  the  right 
hand  side  of  equation  (62)  becomes: 
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i\i2 .  (cos2#  .  cos2  (0-3)  -  1) 

y  .sin#.cos(0-3)  -  iy2  .  cos#  .  (sin2(0-3)  +  sin23) 

-y. sin  (0-3)  -  i\i2  ■  s  i  n2#  .  cos  2  (0  -  3 ) 

Equating  the  second  and  third  elements  of  (63)  and  (65), 
approximations : 


(65) 


and  using  the 


we 


cos(#+A#)  =  cos#/(l  +  tan# .  sinA#) 

cosA0  =  1 

find  approximate  expressions  for: 

sinA0  =  y . tan# . cos (0-3)  -  iy2 . { ( l+2tan2#) . s in2 (0 -3)  +  sin23) 

sinA#  =  -  y.sin(0-3)  -  iy 2 . s in2# . cos2 (0-3) 


(66) 


12.2  Maximum  value  of  bores ight  errors 

For  a  fixed  value  of  y,  and  for  positive  values  of  #,  we  find,  by 
differentiation  with  respect  to  0,  and  subsequent  inspection,  that: 

(a)  s inA#  has  its  maximum  positive  value  when: 


(0-3)  =  -n/2 

3  =  any  value 


(b)  sinA0  has  its  maximum  positive  value  when: 


3  =  -  u/4  or  3tt/4 


sin(0-3)  = 


tan# 


2y  .  ( l+2tan2#) ' 


l  +  2y2 . ( l+2tan2#) 


2  \  i 


tan  # 


(67) 


Figure  5  shows  A9  and  A#  plotted  as  functions  of  y,  for  various  values  of 
#,  and  with  these  values  of  3  and  of  (0-3)-  The  latter  is  also  plotted  in 
the  figure. 

Note  that: 

(1)  3  is  the  angle  of  the  resultant  axis  of  tilt,  due  to  the  combined 
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action  of  ri  and  i p. 

(2)  for  a  fixed  value  of  (r)2+^2),  the  maximum  value  of  AO  occurs  when  ri 
and  < p  are  equal  in  magnitude;  that  is  &  =  -u/4  or  3ir/4. 

(3)  for  a  fixed  value  of  (n2+<J>2),  the  maximum  value  of  A </>  is  (ti2+^2)^, 
and  is  independent  of  $  or 

(4)  the  maximum  value  of  A$  occurs  when  0  =  (5-tt/2;  that  is,  the 
bores  ight  is  normal  to  the  effective  tilt  axis. 

(5)  the  maximum  value  of  A9  occurs  when  0  is  in  the  range 

(3  to  (f5-ir/4)),  depending  on  the  values  of  y  and  <t>  (figure  5(b)). 

(6)  if  the  tilt  of  the  azimuth  axis  is  due  to  a  single  spiral 

deflection,  then  y  =  A  =  n  = 

max  max 

(7)  if  the  pitch  and  roll  errors  are  uncorrelated  sinusoidal  motions, 

then  4>  and  if  can  occur  simultaneously  and  the  value  of  y  is 
max  max  3 

X  XX 

increased  by  22  ;  that  is  y  =  2 2  . ri  =  2  .<J> 

max  max 

(8)  A0  is  not  zero  when  <t>  is  zero,  due  to  second  order  effects. 


13.  EFFECT  OF  TILTED  ELEVATION  AXIS 

13.1  General 

We  wish  to  examine  the  first  order  effect  of  incorrect  alignment  of  the 
elevation  axis  on  the  borcsight  errors  and  on  the  motor  axis  torques. 

If  0  ,  the  origin  of  the  e-system  coordinates  fixed  to  the  elevating 

assembly,  is  displaced  from  its  nominal  position,  the  effect  is  equivalent 
to  a  variation  in  the  coordinates,  X1}  Zj  and  0,  and  can  be  evaluated 

accordingly . 

In  the  remainder  of  this  Section,  we  consider  the  misalignment  in  the 

direction  of  the  elevation  axis,  0  Y  . 

e  e 

13.2  Misalignment  matrix 

Suppose  the  elevation  axis,  0^  Y  ,  is  rotated  through  the  small  angle,  D8, 
about  the  normal  through  0^ ,  parallel  to  OZ^;  followed  by  rotation  through 
the  small  angle,  Dn ,  about  the  normal  through  0  ,  parallel  to  the 
X^Y^-plane.  This  is  illustrated  in  figures  1(d)  and  6.  Then,  to  a  first 
order  approximation  in  DO  and  Dti  : 


-  27  - 
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Subsequent  analysis  is  based  on  the  assumption  that  roll  and  pitch 
correction  errors  are  also  small.  In  this  case,  to  a  first  order 
approximat ion : 


A.  B  =  I 


P,  =  P  =  C. (E  +  M.D.P  ) 
b  s  e 


13.3  Boresight  error 

Using  the  same  argument  as  for  Section  12.1,  it  follows  that: 

(C+AC) . (E  +  M. (D+AD) . (r+Ar) .Z)  =  C.(E+D.r.Z)  (70) 

Neglecting  D  ^E/r  and  Ar/r,  as  in  Section  12.1,  and  multiplying  on  the 
left  by  D* 1 . C~  1 : 

(I  +  D'1 ,C_1 .AC.D) . (I  +  D_1 . (M-I) .D) . (I  +  D'1.  AD).Z  =  Z 
Neglecting  second  order  product  terms  between  AC,  (M-I)  and  AD: 

(D_1.C~l.AC.D  +  D_l.(M-I).D  +  D~l.AD).Z  =  0 

By  evaluation: 


s inA0  -  0 


sinAG  -  Dn  ■  tan<t>  -  D0 


13.4  Elevation  torque  changes 

Changes  in  elevation  torque  due  to  the  inclusion  of  the  misalignment 
matrix,  follow  from  equation  (14): 
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AT  =  -  A  ( R  1  .  R  )  .  ■  K  K  i1 

e 


R*  1  .  R 


n  ■  !'  i  A  i  R 


.R)  .P 


J 


m.N(P)  .I)'1  .D.A(R’  1  ;  .  u;  *  1.)  (72) 

/_, 

e 


Key  results  in  the  evaluation  are  given  in  items  (a)  to  (e)  of  Appendix  II, 
and  the  final  results  are  listed  in  the  left  hand  side  of  Table  7.  Note 
that  the  y-component  is  the  change  in  torque  required  of  the  elevation 
servo  motor . 

13.5  Pitch,  roll  and  azimuth  torque  changes 

Using  equation  (16)  to  determine  the  torque  changes  in  the  r-system: 


ATr  =  A(Q) .Te  +  Q.ATe  -  \  m . N (K) . A (Q . R_ 1 ) . G 


m.N(K)  .A(Q.R  VP’) 


(73) 


Key  results  in  this  evaluation  are  given  in  items  (f)  to  (k)  of 
Appendix  II.  Final  results  are  listed  in  the  right  hand  side  of  Table  7. 
Note  that  the  z-component  is  the  change  in  torque  required  of  the  azimuth 
servo  motor. 


F.quat  iori  (20)  is  used  to  derive  A(T  ,,  +  j.T  .  ,  )  from  AT  ,  and  this  is 

1  roll  J  pjtch  r 

list ed  in  T ub 1 o  8 . 

13.b  Discussion  of  results 


From  ('(juat  ion  (71),  both  Dri  and  D0  contribute  a  significant  error  in 
azimuth  bores  ight ,  but  the  contribution  of  DO  is  constant,  and  can  be 
ca 1 ibrated  out . 

From  Tables  8  and  5,  most  changes  in  roll  and  pitch  torques  are  small 
compared  with  components  existing  with  no  axis  tilt.  The  only  changes  not 
so  swamped,  are  the  terms  in  Dn .  (J  -Q  )  in  the  first  two  columns  of 

G  G 

Table  8.  These  may  be  of  significant  magnitude  in  some  applications. 

A  similar  comparison  between  Tables  7,  3  and  2,  shows  that  none  of  the 

changes  in  azimuth  and  elevation  torque  are  swamped  by  non-perturbed 
torques.  The  most  significant  terms  are  those  in  coupled  with  8  and 

V. 
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Thus  On  is  more  significant  in  its  effects  than  DO,  and  requires  closer 
cont  ro 1 . 


14.  EFFECT  OF  NON- INTERSECTION  OF  AZIMUTH,  PITCH  AND  ROLL  AXES 
14.1  Misalignment  matrix 

We  consider  the  case  of  translation  (without  rotation)  of  the  pitch 
correction,  and  azimuth  axes  so  that  they  do  not  intersect  each  other,  or 
the  roll  correction  axis. 

Figure  7(a)  shows  the  pitch  correction  axis,  0  Y  ,  in  a  plane  normal  to  the 

8  8 

roll  correction  axis,  0,  X.  .  0  0,  is  normal  to  both  axes,  and  is  of  length, 

b  b  g  b  &  ’ 

Z2 .  The  coordinates  of  a  single  point,  P,  are  related  by: 


Pb  =  A.(Pg  +  [0  0  Z2  ]  ) 

Figure  7(b)  shows  the  azimuth  axis,  in  a  plane  normal  to  the  pitch 

correction  axis,  0  Y  .  0  0  is  normal  to  both  axes,  is  of  length,  X,,  and 
8  8  s  g  &  .  z . 

is  collinear  with  the  axis,  0  X  .  0  '0  is  of  length,  Y, .  Then: 

s  s  g  g  &  »  2 

Pg  =  B.(Pg  +  [X2  Y2  0]T) 

Combining  these  relations  with  the  approximation  of  equation  (13): 


Pfa  =  (I  +  X) .  (P  +  E2)  -  Y.E2 


E2  =  {X2  Y2  z2]t 


0  0  Ip 


Y  =  0  0  0 


0  0  0 


14.2  Boresight  error 

Using  the  method  of  Sub-section  12.1: 


A.B. (C+AC) . (E  +  (D+AD). (r+Ar) .Z)  +  (A.B-Y).E2  =  A.B.C.(E  +  D.r.Z.)  (75) 


After  multiplication  by  B  l.A  l/r,  and  neglecting  terms  of  second  order  of 
sma 1 Iness : 
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AC.D.Z  +  C.AD.Z  =  — ( I  - Y ) .  E  2  /  r 


(76) 


Thus  for  small  values  of  the  elements  of  E2  (compared  with  r,  the  distance 
to  the  target),  the  boresight  error  is  negligible. 

14.3  Torque  changes  in  the  s-system 


Using  relations  (10)  and  (77)  to  determine  the  change  in  T^  due  to  a  non¬ 
zero  value  of  E2: 


P  =  P  C.P 

s  r 

P'  =  Pb  =  (I+X).P  +  (I+X-Y).E2 

R  =  I  +  X 

L  =  (I+X-Y).E2 


Clearly,  all  terms  involve  T)  or  ip ,  and  none  involve  time-derivatives  of  0. 
Key,  and  final,  results  are  given  in  items  (a)  to  (c)  of  Appendix  III.  The 
z-componont  of  ATg  is  the  change  in  torque  required  of  the  azimuth  motor. 

14.4  Torque  changes  about  the  elevation  axis 

Using  the  relations  (12)  and  (79)  to  determine  the  change 
component  of  due  to  the  non-zero  value  of  E2 : 


in  the  y- 
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P  =  P  =  C. (E  +  D.P  ) 
s  e 

P  =  d'1 .c'1 .P  -  d'1 .e 

e  s 

Q  =  D_  1  .  c"  1 

K  =  -  D'1 .E 

P',R  from  equation  (77) 


>  (79) 


and : 


AT 

e 


Q.  AT 

s 


L 


m.N(K) .Q.R 


‘  dt 


i>  (  ( I+X-Y)  .  E2  } 


S  D"l.C_l.AT  +  M  .D_1.N(E).c'1.(I-X).(X-Y).E2 
s  e 

-  D~ 1 .C_  1 . (AT  +  M  .N(C.E) . (X-Y) .Ea)  (80) 

s  e 


Again,  all  terms  involve  n  or  < \> ,  and  none  involve  time-derivatives  of  <t> . 
The  y-component  of  AT^  is  the  change  in  torque  required  of  the  elevation 

motor.  Key,  and  final,  results  in  its  evaluation  are  given  in  items  (d)  to 
(g)  of  Appendix  III. 

14.5  Torque  changes  in  roll  and  pitch 
Using  relations  (12)  and  (81): 


Qb  =  A .  B  =  (I+X) 

Kb  =  A. [0  0  Z2]T  +  A.B. [X2  Y2  0]T 

=  A.B. [(Xa-f.Zj)  Y2  Z2)T 
P'  =  Pb  =  (I+X).Ps+Kb 


>  (81) 


then  : 
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=  ( I+X) .AT 


m .  N  (K,  )  . 
/  b 


G  +  Kb+  ^|(I+X)'Ps 


AT  -  N(Kl).M.G  -  N(F.2).  }  M.Kl  +  m.(2X.P  +X .  P  ) 

s  b  ;  b  /  s  s 


-  ( I+X) . N{ ( I -X) . K,} .  m.P 
b  s 


Similarly,  using  relations  (12)  and  (83): 


Qg  =  D 


K  =  B.|Xi  V,  0] 


P'  =  Pu  =  (I+X)  .  P  +  K, 

b  s  b 


then : 


■•N<y  G  +  Kb +  5?(<I*x)'ps,xj 


AT  -  N (K  ).M.G  -  N(B~  1 . K  ).  M . K  +  )  m.(2X.P  +X . P  ) 

s  g  8  b  /  _  s  s 


(I+X) . N{ (I-X) .K  }.  \  m.P 
8  / ,  s 


The  x-component  of  AT^ ,  and  the  y-component  of  are  the  roll  and  pitch 

axis  torque  variations.  Key  results  of  the  evaluation  of  equations  (82) 
and  (84)  are  given  in  items  (h)  to  (q)  of  Appendix  III.  Table  9  lists  the 
results  (as  well  as  those  for  azimuth  and  elevation,  derived  previously) 
except  that  terms  in  n  and  4>  have  been  ignored  wherever  they  are  swamped  by 
larger  terms  in  the  same  cell. 

14.6  Discussion  of  torque  changes 


From  Table  1  and  relations  (22)  and  (25)  we  see  that  H/M,  V/M  and  He/M  are 
small  distances  in  most  applications.  Also,  from  equation  (31): 
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I  ti/O2  |  =  \<p/02\  -  (w2/02).n 

max  max  max 


I  il/O  |  =  |  0/0  j 

max  max 


(w/0)  .  T| 


max 


are  small  quantities  in  most  applications.  Hence  the  terms  in  Table  9  can 
be  divided  into  groups  of  first,  second  and  third  order  of  smallness. 


Only  the  terms  in  g  are  of  first  order  and  these  will  probably  determine 
the  maximum  values  of  X2  and  Y2  ■  Those  terms  independent  of  ti  and  </>, 

together  with  the  elevation  terms  in  M  ,  comprise  the  second  order  group. 

These  may  be  significant.  The  remaining  (third  order)  terms  are  most 
likely  negligible  for  small  values  of  T|  . 


15.  NON -NORMALITY  OF  PITCH  AND  ROLL  AXES 


15.1  Misalignment  matrix 


Consider  the  Cartesian  coordinate  system,  suffixed  g',  derived  by  rotation 
of  the  g-system  (fixed  to  the  gimbal  ring)  through  the  small  angle,  D0 , 
about  the  common  z-axis.  Then: 


P  =  M.P  ' 
g  g 


M 


1  -DO  0 

no  i  o 

0  0  1 


\ 


>  (85) 


Suppose  the  pitch-correction  axis  is  now  transferred  from  the  y-axis  of  the 
g-system,  to  the.  y-axis  of  the  g' -system.  The  pitch  and  roll  axes  are 
still  coplanar,  and  intersect  at  a  common  point,  but  this  intersection 
departs  from  normal  by  the  small  angle,  DO.  Clearly  now: 


PL  =  A  .  M .  B  .  P 
b  s 

and  A.B  is  replaced  by  A.M.B  wherever  it  occurs. 

By  evaluation,  neglecting  terms  of  second  order  of  smallness: 


A.M.B  = 


1  -DO  0 
DO  1  n 
n  l 


(86) 


(87) 


Similarly,  by  evaluation,  we  find  that: 
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A.M.B.C  =  (A.B.C)’ 


(88) 


where  the  slash  denotes  replacement  of  0  by  (0+DG)  wherever  it  occurs. 

15.2  Borosight  error 

Following  the  argument  discussed  in  Sub-section  12.1,  it  is  clear  that: 

A.M.B. (C+AC) . (D+AD) .Z  =  A.B.C.D.Z 


i  f  : 


AG  =  -D0 

A$  -  0 


(89) 


Because  AG  is  constant  it  can  be  calibrated  out. 

15.3  Torque  changes 

It  follows  from  relations  (15),  (17)  and  (18)  that  replacement  of  (A.B.C) 

by  (A.B.C)'  ( ic  replacement  of  0  by  (G+DG)  wherever  it  occurs)  accounts  for 

the  changes  in  T  ,  T  ,  T  .  and  T 

e  r  elev  az 

The  changes  in  roll  and  pitch  axis  torques  follow  from: 


1  ro  1  1 


T 


0  0  1  .  M .  B  .  C .  T  ' 
r 


‘roll 


(90) 


T  .  ,  =  T  =  f  0  1  0 1  .  C .  T  ' 

p l tch  sy  r 


[0  1  0]  . (I  +  ( M " 1  - 1 ) ) .M.C.T  ' 


[’  .  ,  '  -  DO.fl  0  OI.H.B.C.T  ' 

p  i  t<.n  r 


'pitch  '  D0/1roll 


(91) 


15.4  1)  i  sciiss  iori 

Thus  the  effect  of  the  non -norma  1 i ty  is  accounted  for  by  replacement  of  0 
by  fQHlG)  in  all  bores ight  error  and  motor  torque  relations.  In  addition, 
the  pitch  axis  receives  a  component  of  the  roll  axis  torque,  proportional 
to  the  angular  deviation  from  normality.  This  crosscoupling  of  torque  is 
the  only  significant  effect  of  the  non -norma  1 i ty ,  and  cannot  be  eliminated 
by  calibration  of  0. 


■  m 
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TABLE  2.  CONTRIBUTION  OF  ELEVATING  ASSEMBLY  TO  T 


•  .1 
•  •  1 
•  u 


»,  '  ,  T 


v  ■  i.'» ■: » ',■« i i’i '.i  w '5 "M1, ' M s. ' VIH' »■' T 1  j'  j 


»'  .  > 

.  • 

k  \ 

•V* 

L  ^  . 

v  V. 
L" 
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Symbo 1 

Sect  ion 

Do  f ined 

Remarks 

« 

r  .  • 
W' 

First  Used 

*  ■ 

a 

4 

T1 

angle  associated  with  H 

—  ^ 

a 

6 

E3 1 

log.  increment  of  spiral  motion 

y ; 

& 

4 

T1 

angle  associated  with  K 

B 

e 

12 

E64 

angle  associated  with  p 

y 

4 

T1 

angle  associated  with  L 

•' 

n 

1 

- 

roll  error  angle 

•  '  *, 

■ 

r 

0 

1 

- 

azimuth  rotation  angle 

1  -J 

i 

•  ■ 

1 

p 

12 

E64 

amplitude  of  spiral  motion 

<t> 

1 

- 

elevation  angle 

</> 

1 

- 

pitch  error  angle 

^  ■  4 

• 

u 

6 

E31 

angular  frequency  of  spiral  motion 

Wo 

6 

E29 

natural  frequency  of  roll  and  pitch 

m 

i- 

• 

servo 

i 

■J. 

4 

“ 

complex  conjugate  operator 

$ 

[> 

1 

3 

indicates  inertial  coordinate  system 

•  v  -  * 

4 

- 

indicates  replacement  of  0  by  (047r) 

L"* '  ■ 

• 

► 

1 

15 

indicates  replacement  of  0  by  (0+D0) 

: 

tt 

3 

- 

indicates  transformed  coordinate  system 

4 


: 


'\j 


i 
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Symbo 1 

Sect  ion 
First  Used 

De  f ined 

Remarks 

c0 

6 

- 

servo  damping  factor 

d 

6 

E29 

ratio  of  inertias 

e 

2 

- 

suffix  for  e-coordinate-system 

f 

6 

E27 

damping  unbalance  factor 

f 

8 

- 

suffix  denoting  forcing 

f 

8 

E42 

servo  unbalance  function 

g 

2 

- 

suffix  for  g-coordinate-system 

8 

3 

- 

acceleration  of  gravity 

h 

4 

T1 

suffix  for  reference  to  H 

k 

4 

T1 

suffix  for  reference  to  K 

1 

4 

T1 

suffix  for  reference  to  L 

m 

3 

- 

mass  of  point  quantity 

r 

2 

- 

suffix  for  r-coordinate-system 

r 

12 

- 

range  to  target 

s 

2 

- 

suffix  for  s-coordinate-system 

s 

6 

E31 

complex  variable 

t 

6 

- 

t  ime 

u 

6 

E27 

servo  stiffness  unbalance  factor 

X 

3 

- 

Cartesian  coordinate 

y 

3 

- 

Cartesian  coordinate 

z 

3 

- 

Cartesian  coordinate 

A 

12 

- 

operator  specifying  a  change  in  argument 

AO 

12 

- 

azimuth  error 

A  0 

12 

- 

elevation  error 

E 

e 

4 

- 

summation  operator  over  elevating 

assembly 

Z 

r 

4 

- 

summation  over  rotating  assembly 

Z 

r  ,e 

4 

- 

summation  over  rotating  and  all 

elevating  assemblies 
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Section  Defined  Remarks 

First  Used 


2 

2 

3 

4 
3 

3 
2 

9 

8 

8 

4 
4 
6 
4 
2 
4 
2 

14 

2 

14 

14 

2 

12 

2 

14 

2 


origin  of  Cartesian  axes 
point  coordinate  vector 
Ell  transformation  matrix 

T1  mass  distribution  parameter 

E5  transformation  matrix 

torque 

superfix  denoting  matrix 
transposition 

torque  about  azimuth  axis 

torque  about  elevation  axis 

motor  torque 

torque  about  pitch  axis 

torque  about  roll  axis 

E27  servo  stiffness  factor 

T1  mass  distribution  parameter 

Cartesian  axis 

E13  pitch-roll  matrix  approximation 

elevation  axis  offset  component 
pitch-roll  axis  offset 
Cartesian  axis 

E74  matrix  approximation 

pitch-roll  axis  offset 
Cartesian  axis 
[1  0  0]T 

elevation  axis  offset  component 
pitch-roll  axis  offset 
suffix  for  b-coordinate-system 
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NOTATION 


Symbo 1 

Sect  ion 
First  Used 

Defined 

Remarks 

A 

2 

El 

roll  correction  matrix 

A 

6 

E31 

complex  amplitude  of  spiral  deflect 

B 

2 

El 

pitch  correction  matrix 

B 

6 

E32 

complex  amplitude  of  spiral  deflect 

C 

2 

El 

azimuth  rotation  matrix 

D 

2 

El 

elevation  matrix 

Dn 

13 

- 

elevation  axis  deflection  angle 

Dn 

16 

- 

azimuth  axis  deflection  angle 

DO 

13 

- 

elevation  axis  deflection  angle 

DO 

15 

- 

deflection  angle  of  pitch  axis 

E 

2 

El 

elevation  offset  vector 

F.a 

14 

E74 

pitch-roll  offset  vector 

F 

3 

- 

point  mass  force  vector 

F 

6 

E27 

servo  damping  factor 

G 

3 

E6 

gravitation  vector 

G 

6 

E32 

servo  transfer  function 

H 

4 

T1 

mass  distribution  parameter 

I 

4 

- 

unit  matrix 

J 

4 

T1 

mass  distribution  parameter 

K 

3 

Ell 

offset  vector  for  transformation 

K 

4 

T1 

mass  distribution  parameter 

L 

3 

E5 

offset  vector 

L 

4 

T1 

mass  distribution  parameter 

M 

4 

T1 

mass  of  structure 

M 

13 

F.68 

elevation  axis  deflection  matrix 

M 

15 

E85 

pitch -roll  deflection  matrix 

M 

16 

E92 

pitch-azimuth  deflection  matrix 

N 

3 

E3 

torque  generation  matrix 
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to  back  about  the  axis  of  azimuth  rotation,  has  been  examined.  Clear  rules 
are  given  for  determining  whether  the  torques  about  the  azimuth,  pitch  and 
roll  axes,  due  to  mass  unbalance  parameters  of  the  nominally  identical 
structures,  cancel  each  other,  or  are  additive.  For  many  of  the  parameters, 
cancellation  occurs  only  if  the  two  assemblies  are  elevated  identically  and 
s imultaneous ly . 

The  results  of  the  analysis,  presented  in  tabular  and  graphical  form,  provide 
a  firm  basis  for  the  design  of  the  servo  drives  and  of  the  procedure  and 
facilities  for  achieving  the  necessary  degree  of  mass  balance. 

The  motion  of  the  stabilised  platform  can  be  described  in  terms  of  a  sequence 
of  small  amplitude  spiral  nutations  of  the  azimuth  axis  about  the  true 
vertical  direction. 

Two  of  the  main  components  consist  of  the  transient  decay  of  an  initial 
deflection  of  the  azimuth  axis  from  vertical.  One  of  these  is  in  the  forward 
direction  (the  same  direction  as  the  azimuth  scanning  rotation)  and  the  other 
is  in  the  reverse  direction.  They  are  mirror  images,  although  executed  at 
widely  different  rates.  The  actual  angular  velocity  of  the  nutations  depends 
on  the  stiffness  and  damping  factor  of  the  pitch  and  roll  servo  loops,  the 
rate  of  azimuth  rotation,  and  the  distribution  of  mass  of  the  combined 
rotating  and  elevating  assemblies. 

The  third  main  component  is  a  non-decaying  (circular)  nutation  in  the  forward 
direction,  and  at  the  same  rate  as  the  azimuth  scanning  rotation.  Its 
amplitude  depends  on  the  degree  of  unbalance  in  the  distribution  of  mass,  as 
well  as  on  the  servo  stiffness  and  damping,  and  the  azimuth  rate. 

Each  of  the  three  main  components  spawns  two  minor  components  proportional  to 
its  amplitude.  One,  at  the  same  rate  and  opposite  direction  to  the  parent 
nutation,  is  proportional  to  the  unbalance  in  the  stiffness  and  damping  of  the 
pitch  and  roll  servo  loops.  The  other,  at  a  rate  and  direction  depending  upon 
both  that  of  the  parent  nutation  and  of  the  scanning  rotation,  is  proportional 
to  another  measure  of  the  unbalance  in  the  mass  distribution. 

Each  of  the  minor  components,  in  turn,  spawns  two  more  minor  components,  ad 
infinitum,  but  the  amplitudes  of  all  these  are  negligible  for  modest  amounts 
of  unbalance. 

Many  imperfections  of  mass  balance  and  axis  geometry  have  an  effect  on  the 
accuracy  of  the  boresight  direction,  but  most  are  negligibly  small,  or  are 
cancelled  out  during  the  calibration  process.  Those  that  have  a  significant 
residual  effect  (and  therefore  require  appropriate  control)  are: 

(a)  the  instantaneous  amplitude  of  the  total  nutation  of  the  azimuth  axis 
about  the  vertical 

(b)  rotation  of  the  elevation  axis  from  its  nominal  position,  in  the 
vertical  plane 

(c)  non-normality  of  the  pitch  and  azimuth  axes. 
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9 

m 


Lk 


T  .  .  =  T  =  [10  0] .B.M.C.T  ' 

roll  gx  1  r 

=  [100] .A. B.M.C.T  ’ 

r 


=  T 


roll 


(95) 


T 

P 


i  tch 


T 

sy 


[0  1  0] ,M.C.Tr' 


=  [0  1  -Dn] . C .T  ' 


=  T 


pitch 


-  Dr)  .T 


az 


(96) 


16.4  Discussion 

Thus  the  effect  of  the  non-normality  is  accounted  for  by  replacing  ri  by 
(ti+Dti )  in  all  boresight  error  and  motor  torque  relations.  In  addition,  the 
pitch  axis  receives  a  component  of  the  azimuth  axis  torque  proportional  to 
the  angular  deviation  from  normality.  This  cross  coupling  of  torque  is  the 
only  significant  effect  of  the  non-normality,  as  it  cannot  be  eliminated  by 
the  calibration  of  n . 


17.  CONCLUSION 

The  four-axis  stabilised  director  system,  while  mounted  in  the  normal  erect 
position  on  a  firm  foundation,  has  been  analysed  to  determine  the  driving 
torques,  boresight  errors,  and  motion  of  the  stabilised  platform  caused  by  a 
smal 1 : 

(a)  deflection  of  the  stabilised  platform  from  the  horizontal  position 

(b)  unbalanced  mass  distribution  on  the  rotating  and  elevating  assemblies 

(c)  unbalance  between  the  gain  parameters  of  the  roll  and  pitch  correction 
servo  loops 

(d)  deviations  from  the  ideal  geometrical  relationships  between  the  axes. 

Techniques  are  presented  for  factorising  the  complex  matrix  relationships  in 
such  a  way  as  to  greatly  simplify  their  evaluation. 

The  approximation  of  trigonometric  functions  of  the  roll  and  pitch  angles  by 
first  order  functions  of  the  angles  themselves  (justified  by  the  assumption 
that  these  angles  are  always  small)  greatly  simplifies  the  interpretation  of 
the  results  of  the  analysis  by  reducing  significantly  the  number  of  terms 
contributing  to  each  relationship.  Conclusions  drawn  in  respect  of  a 
particular  application  should  be  carefully  checked  for  validity;  if  necessary 
using  the  accurate  relationships.  This  is  particularly  important  if  the 
analysis  were  to  be  extended  to  the  case  of  the  director  system  mounted  on  a 
pitching  and  rolling  vehicle. 

The  case  of  two  separately  elevatable  assemblies,  mounted  symmetrically  back 
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16.  NON-NORMALITY  OF  PITCH  AND  AZIMUTH  AXES 


16.1  Misalignment  matrix 


Consider  the  Cartesian  coordinate  system,  suffixed  s',  derived  by  rotation 
of  the  s-system  (fixed  to  the  stabilised  platform)  through  the  small  angle, 
Dn  ,  about  the  common  x-axis.  Then: 


P  =  M.P  ' 
s  s 


M 


1  0  0 
o  i  -Dn 
o  Dn  l 


\ 


>  (92) 


/ 


Suppose  the  azimuth  axis  is  now  transferred  from  the  z-axis  of  the  s-system 
to  the  z-axis  of  the  s' -system.  Then  the  azimuth  and  pitch  axes  intersect 
at  an  angle  departing  from  normality  by  the  small  angle,  Du.  Clearly  now: 


P,  =  A  .  B  .  M .  C  .  P 
b  r 


(93) 


By  evaluation,  neglecting  terms  of  second  order  of  smallness: 


A .  B .  M 


1  0  ip 

o  1  -(n+Dn) 

ip  (n+Dn)  1 


\ 


>  (94) 


=  (A. B) ' 


where  the  slash  denotes  replacement  of  ri  by  (n+Dn)  wherever  it  occurs. 

16.2  Boresight  error 

From  equation  (94),  the  effect  of  Dn  is  clearly  the  same  as  for  the  same 
magnitude  of  n,  as  determined  in  Sub-section  12.1. 

16.3  Torque  changes 

It  follows  from  relations  (15),  (17)  and  (18)  that  replacement  of  (A.B)  by 

(A.B)'  ( ie  replacement  of  n  by  (n+Dn)  wherever  it  occurs)  accounts  for  the 

changes  in  T  ,  T  ,  T  .  and  T  .  The  changes  in  roll  and  pitch  axis 
e  r  elev  az 

torques  follow  from: 


c\ "  *x- 

guis-  t2'  ‘Z- 


-  il  - 
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APPENDIX  I 


KEY  RESULTS  IN  THE  EVALUATION  OF  T  AND  T 

e  r 


Making  the  substitution: 


a  —  sin0 


b  =  cos0 


c  =  sin 0 


d  =  cos0 


(a)  C_l.X.C 


0  -  Ca  -  n -b . </#)  [ 


0  -(b.n+a.0>) 


(a.n-b.0)  (b.n+a.00 


=  i(n+j.0).e  To  0  -j~]  +  Kn-j -0)  .e“i0.  Fo  0  j 


0  0-1 


0  0-1 


j  1  0 


-j  1  0 


Note  that  the  second  term  is  the  complex  conjugate  of  the 
first  term. 


(b)  C_1.X.C  =  i(h+j.0).e‘J  .0.  000  +  Conjugate 

0  0  0 


1  -j  0 


(c)  C'FC  =  0.  0  -1  0 


1  0  0 


0  0  0 


(d)  C*1 .C  =  0.  0  -1  0  +  02  .  -1  0  0 


1  0  0 


0  0  0 


0-10 


0  0  0 
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(e)  D_l.D  = 


0  0-1 

0  0  0 

1  0  0 


(f)  R- 1 . R 


D'1.C'1.(I-X).{X.C.D  +  (I+X).(C.D  +  C . D) } 
D_1 .C'1 .X.C.D  +  D"1.C'*.C.D  +  D_ 1 . D 


=  i(n+j-0)-e 


-j0 


0  c  -  j 
-c  0  -d 
j  d  0 


+  0. 


+  Conjugate 


0 

-d 

0 

+  0- 

~0 

0 

-1 

d 

0 

-c 

0 

0 

0 

0 

c 

0 

1 

0 

0 

(g)  )  m.N(P).R_1.R.P  = 


-<P)  -e 


■j0 


L_ 


(Je+Qe).d  +  K^.cosCtf-fJ^)  +  j.L^.cosy^ 

-1^.  cos  (0+3^)  -  2j.Je 

-(J  +Q  ).c  +  K  .sin(0-B  )  +  j.L  .sinlf 
e  e  e  e  e  e 


+  Conjugate 


(h) 


+  0. 

(J  +Q  ) .c  +  K  .sin(0-&  ) 
e  e  e  e 

+  0. 

L  .  cosJf 
e  e 

-L  .sin(0+Y  ) 
e  e 

-2J 

e 

'  K  cos(0-B  ) 
e  e  e  e 

L  .sin2f 
e  e 

N(P) .D'1 

-c.V  -H  .sina 

e  e  e 

d.V 

e 

H  .sin(0+a  )  0 

e  e 

-H 

e 

cos(0+a  ) 
e 

-d.V 


H  .cosa 
e  e 


-c.V 
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(i)  D . R* 1 . G  =  C'1.(I-X).G 

+  Kn+j  .1 l/)-e  j0. 


0 

0 

g 


(j)  D . R~ 1 .  L 


C'1.(I-X).{X.C.E  +  2X.C.E 
{C_1.X.C  +  2C_1 .X.C  +  C"1 


=  i(Ti+j.i^).e 


■j0 


-j-Zx 

-Zi 

j-Xi 


+  Cn+j - 0) - e  ^  0 . 8 . 


0 

0 

Xi 


+  0. 

0 

+  92  . 

~-xr 

X! 

0 

0 

0 

(k)  N(E) 


0  -Z1  0 

Zx  0  -X! 

0  Xi  0 


(1)  M  .E  +  V  m.D.P 

e  L 

e 


M^.Xi  +  He-cos^^ 


j.g  +  Conjugate 

g 

0 

i  +  (I+X).E> 

•  C}.E 

■  Conjugate 

+  Conjugate 


M  . Zj  +  H  .  sin^ 
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(m)  \  m.D.P  =  #.  -H  .sin #, 

/  e  e  h 


H  .cos #, 
e  h 


(n)  m.D.P  =  #.  -H  .sin#,  +  #2.  -H  .cos#, 

/  e  eh  eh 


H  .cos#, 
e  h 


-H  .sin#, 
e  h 


(o)  N(E)  .C" 1  .X.C  =  i(Ti+j.#).e”J  .  0  0  Zi  +  Conjugate 


-j.Xi  -X,  -j.Zx 
0  0  -Xx 


(p)  N(E).C  1.X.C.  =  i(n+j.#).8.e  ^  .  0  0  0  +  Conjugate 

-Xi  j.Xx  0 
0  0  0 


(q)  N ( E ) . C " 1 . C  =  0.  -Z,  0  0  +  02  .  0  Zx  0 


0  -Zx  0 


X,  0  0 


■Z,  0  0 


0  -X,  0 


fr)  N(E).C  '.X.C  =  i  (n+j .  .  e  ^  .  0  0  Zx  +  Conjugate 

-j.X,  -x i  - j . Z 1 
0  0  -X, 


(s)  N(E ) . C  .C  =  0.  -Zj  0  0 


0  -Z,  0 


X,  0  0 


m . N(K) . Q . R- 1 .G  =  M  . N(E) . c" 1 . ( I -X) . G 
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m . N(K) . Q . R" 1 . P1  = 


N(E) . (C_ 1 . X . C  +  2C_1 .X.C  +  C_1.C).'  M  .E  +  \  m.D.P 

e  /  e 


+  2N(E).  (C_1.X.C  +  C^.CjA  m.D.P  +N(E).\  m.D.P 

’  e  /  e 


R  =  C_1.X.C  +  C_1.C 

m.N(P) . (R- 1 . R . P  +  P)  = 


0. 

-L  .cosy 
r  r 

+  i(Tl+j  •’)')  -e  . 

J  +Q  +K  .  cosfS  +j.K  .  sinft 
r  r  r  r  J  r  Kr 

-L  .sinT 
r  r 

-j(J  +Q  )+j-K  .cosP  -K  .sinft 
r  r  r  r  r  r 

2J 

i  .  L  .  sinfr  -  L  .  cos2T 

r 

r  r  r  r 

+  Conj”t>ate 


m.N(P) 


0 

V 

r 

-H  .sina 
r  r 


-V 

r 

0 


H  .sina 
r  r 

-H  .cosa 
r  r 


H  .cosa  0 


G  = 


C~ 1 . ( I -X) . G 
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APPENDIX  II 

KEY  RESULTS  IN  THE  EVALUATION  OF  EQUATIONS  (72)  AND  (7.3) 
(a)  A(R_ 1 .R)  =  D_1 .M_1 .C_l . (C.M.D  +  C.M.D)  -  D_l .C_l . (C.D  +  C.D) 


=  D* 1 . {C~ 1 . C . (M-I )  -  (M-I).C''.C).D 


=  Dn .  9 .  0  0  1 


0  0  0 


-1  0  0 


(b)  '  m. N(P) . A(R" 1 . R) . P  = 


Dri .  0  .  -L  .cosJT  I 
e  e  | 


L  .  sin2f 
e  e 


(c)  D . A(R~ 1 )  =  D . D~ 1 . M_ 1 . c" 1  -  D.D~ 1 .C" 1 


=  -  (M-I ) . c" 1 


(d)  D .  A(R_  1 ) .  G  =  Dri .  g .  |  0 


(e)  D.A(R-l)-L  =  -  (M-I) .C_1 .C.E 


0.X,.  D0  +  02.X,.  0 
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APPENDIX  III 


(a) 


r,e 


KEY  RESULTS  IN  EVALUATION  OF  SECTION  14 
[H.cosB  H.sinG  V]^ 


(b)  (X-Y).Ej  =  (0  -n.Z2  (r\ .  Y2  -ip  .  X2 )  ] T 


(c) 


AFs  =  ^ ' 


-V.Z2  -  Y2.H.sin0, 


Y2  . H , cos0^ 


Z2  .  H .  cos0^ 


+  <P- 


X2 . H . s in0 


-X2 . H . cos0 


(d)  [0  1  0].D  *.C  1  =  [-sin0  cosG  0] 

(e)  C.E  =  [Xi-cosO  Xj.sinG  ZX]T 


(f)  [0  1  0 ] . D~ 1 . C " 1 . N(C . E )  =  [ Z i . cos0  Z j . s in8  -X2 ] 


(g)  ATelev  =  V.{(V-Me.Z1).Z2.sin0  -  Y2 . (Me .X, -H. cosa) } 

+  ip.X2.(M  .Xj-H-cosa) 
e 

(h)  ’  m.P 

s 

r,e 


H .  cosG, 


H .  s in0, 


(i) 


m.P 


■0.H.sin0  -  0.H  .cosO.sin0. 

he  h 

O.H.COS0,  -  0.11  .  sin0.sin0. 
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r  ,e 


00  Kb 

(1)  [1 

(m)  [  1 

(n)  [1 

(o)  [0 

(P)  [0 


m .  P 


-O.H.sinQ.  -  0j.H.gos0  +  20.0.  It  .  s  inO  .  s  in0, 
h  h  e  h 

-  0.11  . cosO.sin0,  -  0J.H  . cos0.cos0, 
e  he  h 


G.H.cosQ,  -  G2.H.sin0,  -  2G.0.H  . cos0.sin0, 
h  h  e  h 

-  0.H  .sinG.sin0,  -  02.H  .sinG.cos0. 
e  he  h 


L 


0 


x2 

Y2  -  .7. i 

Z2  +  n  .  Y 2  -  ip.Xj 


0  0].N(Kb) 

0  0] .N(E2).X 
0  0] . C I+X) . N{ ( I -X) . Kb> 
1  0].N(K  ) 

1  0].N(B"l.K  ).X 


(0  -(z2+n.Y2-0.x2)  (Yj-ti.Zj)] 


I-*.y2 

-(Z2-*.X2) 

Y2] 

[-*.y2 

-(Z2-*.X2) 

y2] 

[-*.X2 

0  -X2] 

[*.X2  - 

T1  .  X2  "Xj  ] 

(q)  [0  1  0|.(I+X).N{(I-X).K  } 


[0  -T| .  X2  -X2  I 
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